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Abstract 

We have studied how lepton pairs from decays of heavy-flavoured mesons produced 
in pA collisions can be used to determine the modifications of the gluon distribution 
in the nucleus. Since heavy quark production is dominated by the gg channel, the 
ratio of correlated lepton pair cross sections from DD and BE decays in pA and pp 
collisions directly refiects the ratio = fflfg- We have numerically calculated the 
lepton pair cross sections from these decays in pp and pA collisions at SPS, RHIC and 
LHC energies. We find that ratio of the pA to pp cross sections agrees quite well with 
the input R^. Thus, sufficiently accurate measurements could be used to determine the 
nuclear modification of the gluon distribution over a greater range of x and than 
presently available, putting strong constraints on models. 



^ kari.eskola@phys.jyu.fi 
^ vesa.kolhinen@phys.jyu.fi 
^vogt@lbl.gov 



1 Introduction 



Inclusive differential cross sections of hard scatterings in proton-nucleus collisions at 
high energies are, to first approximation, computable assuming factorization. For pro- 
cesses not sensitive to isospin effects, such as heavy quark production, these cross 
sections can be expressed as 

where alQ"^ , Xi, X2)ij^a+x are the perturbatively calculable partonic cross sections for 
producing parton a at scale ^ ^qcd? X2 are the momentum fractions of the 

partons involved in the hard scattering, while ff and //^ are the distributions of parton 
i in a free proton and a nucleus with mass number A respectively. In the high-Q^ limit, 
possible multiple scatterings of individual partons ||1|, 0] are neglected but the parton 
distributions will be modified in the nucleus 0, Equation (|l|) may be generalized to 
nucleus-nucleus collisions but we focus only on pA collisions in this paper. 

The measurements of the structure function in deep inelastic lepton-nucleus 
scattering (DIS) clearly show that, compared to deuterium, D, the ratios Rp^ = 
F-f/F^, differ from unity. As nuclear effects in D are negligible, this indicates that 
the parton distributions of bound nucleons are different from those of a free proton, 
Rf{x,Q'^) = //^(x, (x, Q^) 7^ 1. The nuclear effects observed in F^ are usually 

categorized according to the different regions of Bjorken-x as: shadowing, Rp^ < 1, 
at x<0.1; antishadowing, i?^ ^ 1' 0.1<x<0.3; the EMC effect, i?^ ^ 1' 
0.3<x<0.7; and Fermi motion, Rp^ > 1, at a; — > 1 and beyond. For a review of the 
measurements and some models of nuclear effects, see Ref. j^. 

The scale evolution of the nuclear parton distributions is, to a first approxima- 
tion, similar to that of the free proton. Thus at scales > Ql ^ ^qcd the evolu- 
tion of //^(x, Q^) can be described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations [|] if the initial conditions are given at a fixed initial scale 
Qo- Such a DGLAP analysis of the nuclear parton densities |^, ^, 0| assumes that the 
nuclear effects factorize from the parton densities for Ql^l GeV^ and the subsequent 
Q^-evolution of the modifications follows the DGLAP equations, neglecting ~ 
corrections |||, ^j. The nuclear modifications at Ql, Rf{x, Ql), can then be determined 
based on constraints provided by the measurements of F^ in DIS, the Drell-Yan cross 
sections in pA collisions |0, |lT] and conservation of momentum and baryon number. 
This procedure was carried out in Ref. In the present work, we utilize the EKS98 
parametrization of the nuclear effects Q^), obtained from Ref. 0, released for 

public use in Ref. 0. 

There remain uncertainties, especially in the determination of the nuclear gluon 
distributions. Figure |l| shows the gluon modifications, Rf{x,Q^), for a heavy nucleus, 
A = 208, as obtained from a DGLAP analysis |0, ^. So far there are no direct 
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Figure 1: Scale evolution of the ratio R^{x,Q'^) for an isoscalar nucleus A=208 according 
to EKS98 I, i- The ratios are shown as functions of x at fixed values of equidistant in 
logg^: 2.25 GeV^ (solid), 5.39 GeV^ (dotted), 14.7 GeV^ (dashed), 39.9 GeV^ (dotted- 
dashed), 108 GeV^ (double-dashed), and 10000 GeV^ (dashed). The regions between the 
vertical dashed lines show the dominant values of X2 probed by muon pair production from 
DD at SPS, RHIC and LHC energies. 



constraints on the nuclear gluon distributions. These can, however, be constrained 
indirectly through the scale evolution of F^(x, Q^) [^, |13|. The only measured 
dependence is of the ratio F2^/F2 at 0.02 < a; < 0.2. At other values of x, only 
the requirement of stability of the evolution can be used in the perturbative analysis. 
Finally, a global limit on the nuclear gluon distribution is provided by momentum 
conservation v^^hich effectively determines the level of antishadowing in Rf{x, Ql). More 
measurement-based constraints are thus vital for reducing the uncertainties in the 
nuclear gluon distributions. Heavy quark production is one such process. 

The dominant subprocess of heavy quark production in hadronic or nuclear colli- 
sions is gg QQ. The produced heavy quarks fragment into mesons which may decay 
semi-leptonically. Lepton pairs formed from these decay leptons thus carry direct in- 
formation about the input nuclear gluon distribution. Previous studies can be found 
in Refs. |15|. At that time, however, R'^ was not yet as well constrained as in the 



EKS98 parameterization. 

There are clear advantages for studying the nuclear gluon distribution in pA inter- 
actions. In pA studies, any possible complications of dense collective matter produced 
in ultrarelativistic AA collisions are avoided. Also, the number of uncorrelated pairs 
where e.g. the D and D originate in two independent partonic subcollisions remains 
rather small. 
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The purpose of this paper is to study to what extent lepton pair production from 
heavy flavor meson decays reflects nuclear effects on the gluon distributions in pA 
collisions, especially in the ratio ctP^/ctPP. As the input nuclear gluon modification, we 
will use the EKS98 gluon parameterization shown in Fig. [|. Both correlated lepton 
pairs, where the QQ pair is created in one subcollision, and uncorrelated pairs are 
considered. To probe the ratio Rf{x, Q"^) over as wide an x range as possible, we study 
pA and pp collisions at the SPS (^i = 17.3 GeV), RHIC (v^ = 200 GeV) and LHC 
{^/s = 5500 GeV) assuming the same energy as in AA collisions for both pp and pA 
studies. Lepton pair production from DD decays at the SPS (NA60) will probe the 
antishadowing region while at RHIC (PHENIX, STAR) and LHC (ALICE, CMS) the 
cross sections are predominantly sensitive to the gluon properties in the shadowing 
region, as indicated in Fig. [1|. 

It should be noted that the present work focuses only on the best-case scenario in 
which we assume that the leptons from heavy meson decays can be identified through 
vertex displacement. We make no attempt to simulate the background processes in 



detail here. This study, together with Refs. [IJ, |15|, should provide further motivation 
for performing more detailed event simulations in real detectors. Our hope is that 
the lepton pair cross sections measured in the future could be used to constrain the 
remaining uncertainties in the nuclear gluon distribution in a model-independent way. 



2 Formalism 

2.1 Lepton pair production in pA collisions 

We first present the expression for the differential cross section of lepton pair production 
in pA collisions. To simplify notation, we refer to generic heavy quarks, Q, and heavy- 
fiavoured mesons, H. In general, the lepton pair production cross section is 



,^ , . J d^Pid^Pi J d^pHd^p^SiM^j - M{pi,pj)) 5{yii - y{puPi)) 

II ^11 

rfr^-'+^(pW) rfr^-^+^(%) rfaP^-^^+^ 

d'^Pi d^pi d^pnd^pjj 

X6'(l/mm < yi,yj < yma.x)0{4> mill ^ (plj 4^1 ^ 0max 

)• (2) 

where M^j and are the mass and the rapidity of the lepton pair, defined as 



M{pi,Pi) = ^Jipi+Pi)', (3) 

y[pi,pi) 2'''[{E, + Ej)~{p,, + pj^)) ■ 
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The decay rate, dT^^'-^^ (pn) /d^Pi, is the probabihty that meson H with momentum 
Ph decays to a lepton / with momentum p/. The 9 functions define single lepton rapidity 
and azimuthal angle cuts used to simulate detector acceptances. 

Using a fragmentation function Dq to describe quark fragmentation to mesons, the 
HH production cross section can be written as 

X ^'''(PH-2ijf<3)<5'"fe--2R3) (6) 

where z is the fraction of the parent quark momentum carried by the final-state meson. 

Quark fragmentation has been described by several different models. A fragmen- 
tation function smears the quark momentum into the momentum of the final-state 
meson during hadronization. Typical fragmentation functions are a delta function or 



the Peterson function ||T^]. In this study we have used the delta function but have also 
checked our results using the Peterson function. Delta function fragmentation simply 
assumes that the quark and meson momenta are identical. 



D^{z) = 6{l-z) . (6) 

On the other hand, the Peterson function assumes that the meson momentum is 
smeared according to 



N 

z{i-{i/z)~Q/ii-z)y 



DS (^) = .... ..... , ... • (7) 



The normalization, N, is fixed by the requirement J2h Io '^^ Dq{z) = 1. The peak of 
the fragmentation function is at z ~ 1 — 2eQ with a width eg ~ {mq/mQ)'^ We 
have assumed e.g. that all the c quarks fragment into D mesons so that if = D for 
charm. 

The hadronic heavy quark production cross section per nucleon in pA collisions can 
be factorized into the general form 

The parton distributions /f and ff, described in Eq. (|l|), are evaluated at the scale 
Q"^ ~ m"^^. To lowest order, LO, the partonic cross section is 

^ ^ da'^^^'^ s da'i-^^'^ 
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SPS (A = 208) 


RHIC {A = 197) [/ib] 


LHC {A = 208) [/ib] 


PP 


1.188 


143.7 


6137 


~cc 
^pA 


1.282 


140.1 


4826 


~hl 
PP 


1.674 X 10-5 


1.412 


251.9 




1.464 X 10-5 


1.513 


226.7 



Table 1: The total QQ production cross sections per nucleon in pp and pA interactions. 

where pi and p2 are the four momenta of the incoming partons. The lowest order 
convolution of the parton distributions with the subprocess cross section is written as 



E [/.^(^i,Q')/#(^2,Q^) + /|(xi,g^)/,^(x2,Q^)]|^,,^QQl ) . (10) 

q=u,d,s 



The squared matrix elements \M.„g^QQ\ and | A^^^^qqI can be found in e.g. Ref. [|T7 

Integration of Eqs. (|) and (|5j) over the total phase space gives the normalization 
^pA^a ^ Q2^pA^HH ^ Q2^pA^QQ^ ^^^^^ q ^ Yi/T. The total cross section for lepton 
pair production is then equal to the total cross section of heavy quark pair production 
multiplied by the square of the lepton branching ratio. The following expressions for 
the double differential distribution and the total cross section of QQ production per 
nucleon will be useful for our later discussion: 



1 d^pA^QQ+X. 



da^^ _ 1 1 

da^^ 



-EqEq- 



(11) 



= J dpldy,dy2j:xjr{xuQ')x2ff{x2,Q')^^ , (12) 

where Xi^2 = rriT / y/s {e'^^'^ + e^^^) and m'^ = Pt + "^q- The rapidities of the heavy 
quarks are 1/12. The total QQ cross sections per nucleon in pp and pA interactions are 
presented in Table |l]. 

2.2 Correlated and uncorrelated pairs in pA collisions 

In pA collisions, a number of uncorrelated lepton pairs are also produced along with the 
correlated lepton pairs. Since, in most cases, one would prefer to study either correlated 
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or uncorrelated pairs alone, these events must be separated. One method of separation 
is like-sign subtraction because like-sign pairs are all uncorrelated: /+'s originate from 
e.g. c quarks whereas /~'s come from c quarks. Another possible method is to use 
azimuthal angle cuts. At lowest order, correlated QQ pairs are produced back-to- 
back. Thus, leptons produced from these pairs tend to have a larger relative azimuthal 
angle than uncorrelated ones. With a suitable azimuthal cut the relative amount of 
uncorrelated pairs could possibly be reduced. However, the focus of our study is not 
how to separate the pairs but an estimate of their rates. 

In order to calculate the number of correlated and uncorrelated pairs, we assume 
that all subprocess interactions are independent. In this case, the number of interac- 
tions can be described by a Poisson distribution. In our determination of the relative 
number of correlated and uncorrelated pairs, it is convenient to consider the impact- 
parameter integrated QQ cross section, 

■ £h (1 - e-^-(^)'^'='^) = fd'h± ^QQ(b)^^''"^^'^ , (13) 

N=l ^■ 

where b is the impact parameter. The number of QQ pairs produced as a function 
of impact parameter is A^QQ(b) = TA{h)a^^ with T4(b) = / dzpAih,z) the thickness 
function of a nucleus with density pA(b, z). 

Multiplying the right hand side of Eq. (0) by the unit normalized QQ cross section, 

1^ = f J= / d%d^i%-^^^Y , (14) 

and differentiating with respect to d^pgd^pq we find 

N 



^pA^QQ+X 



N 

X ( E ^^'Hpq - M^^'Hpq - PQk)) ■ (15) 
j,k=i 



d^Pqd^PQ 



We have written the power A^ in Eq. (|T4|) as a product in order to separate the momenta 
Pi of different subcollisions. We have also written the delta functions as a sum over 
j and k in order to indicate to which subcollision they belong. When considering 
correlated pairs, both quarks come from the same subcollision so that j = k. In the 
case of uncorrelated pairs, j ^ k, the Q and Q are picked from different subcollisions. 
Thus for correlated pairs, j = k, we have 

^ pA^QQ+x oo^ N^j^(h)e-^QQ^^^ 1 



rf^bE 



d^Pqd^PQ J N\ a^Q d^pqd^p^ 

^ da^^ 
d^Pgd^Pg 
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Similarly, for uncorrelated pairs, j ^ k, 



d^PQd^PQ J ki N\ ^ ^ aQQ d^pg aQQ d^p^ 

We use the definition of the nuclear overlap function, T^bO^) = J dsTA{s)TB(h — s), 
with A = B to obtain the integral over b in Eq. (P), / d'^hTjih) = Taa{0) = 
A'^/{7iR\). The last equality is obtained by assuming diffuse Woods-Saxon distribu- 
tions for the nuclear density. Thus, for correlated pairs we find 

J pA^QQ+X J~QQ 
""corr ^ j{_ZZ__ 

d^Pqd^pQ d^pQd^Pg ' 

while for uncorrelated pairs 

rfaPt7««+^ A^ da'^^da'^'^ 



d^Pqd^Pg TiR\ d^pQ d^pg 



(19) 



The double differential cross section for the pairs, the right-hand side of Eq. ([Tsp, was 
defined in Eq. ([Tl|) . The single Q cross section, da^^ / d^pq^ is the integral of Eq. (|Tl]) 
over d^pq. Dividing Eqs. ([T8|) and (|19|) by A gives the QQ production cross section 
per nucleon. 

The above formalism can be extented to AB collisions by simply replacing the 
thickness function Ta in Eq. (|13|) by the nuclear overlap function Tab- 

As seen in Eq. (|T^), the uncorrelated pair cross section is proportional to the product 
of two independent single heavy quark cross sections. The heavy quark cross sections 
are, in turn, proportional to the lepton cross sections. This proportionality to the lepton 
cross sections could be utilized to study the ratio of cross sections of uncorrelated lepton 
pairs from pA and pp collisions. This uncorrelated ratio would probe the square of the 
nuclear modifications of gluon distribution function, {RfY, instead of Rf. Because the 
uncorrelated pair cross section is the product of independent single lepton distributions, 
lepton pairs produced by leptons from two different nucleon-nucleon collisions can be 
used to simulate uncorrelated pairs. 



3 Numerical inputs 

We calculate the lepton pair cross section in Eq. (||) by generating lepton pairs in 
proton-nucleus collisions by Monte Carlo and binning them according to their mass 
and rapidity. We calculate (xi), {X2) and {Q) in each {yij,Mfi) bin. The {X2) regions 
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probed are indicated in Fig. |T]. For simplicity, since isospin effects do not play a role 
in gg QQ, we have assumed that the nuclei consist of A protons. We have used the 
MRST LO (central-gluon) parton distribution functions ||T8| from PDFLIB [|T^ along 
with the EKS98 [|| nuclear modifications. The quark masses are rric = 1.2 GeV and 
fUb = 4.75 GeV while the meson masses are = 1.8693 GeV and rriB = 5.278 GeV. In 
cc production we have chosen the factorization scale to be twice the quark transverse 
mass, = 4mf.^, because the lowest values of m^^ might otherwise be below the 
minimum Q'^ of the parton distribution set. Since the transverse mass is larger in bb 
production, we use = m^^. 

The c quark fragments into D~^, and Dg mesons with respective relative weights 
1:1:0.3 Here we neglect Dg mesons and allow c quarks to only fragment to or to 
Z)° with equal probabilities which should not significantly affect the results. Similarly: 
c-^ D- and b ^ B+ and 5^; and b ^ and . _ 

In order to obtain sufficiently high statistics, we generated 2 x 10^ QQ pairs and 
let each pair fragment 10^ times. We used JETSET/PYTHIA subroutines [|^] for 



meson decays to leptons. To further improve statistics and reduce the runtime of the 
calculation, we allowed leptonic decays only. The final normalization is obtained using 
the leptonic branching ratios. The total QQ production cross sections per nucleon 
needed for normalization are given in Table 1. 

We also generated uncorrelated lepton pairs from two QQ pairs which then sepa- 
rately fragment and decay. The lepton pairs originating from different quark pairs are 
the uncorrelated ones. This is the origin of the factor of two in the previous paragraph. 

We then binned all correlated and uncorrelated e^e^ ^ e^/i''', and n^^i^ pairs. We 
accept only leptons which originate directly from DD and BB meson decays and neglect 
the background from secondary or tertiary leptons from e.g. — > K^X — > fi^X'. 

Our pp and pA results are calculated for SPS {^/s =17.3 GeV), RHIC (v^ =200 
GeV) and LHC (^5=5500 GeV), using the same center of mass energies as in AA colli- 
sions to probe the same x values in all collisions. At the SPS and LHC, we assume a Pb 
(^=208) target and an Au (A=197) target for RHIC. The new NA60 ||2l| experiment 
at the SPS makes it possible to detect muon pairs in the range < 77^ < 1. Since 
the low energy makes bb production negligible, only cc production is studied at the 
SPS. At RHIC, the PHENIX detector |^] is designed to measure ee, e/i and /x/x pairs. 
The detector has two electron and two muon arms. The acceptance of electron arms 
is —0.35 < rje < 0.35 in pseudorapidity and ±(22.5° < 0e < 112.5°) in azimuth. The 
forward muon arm covers the pseudorapidity region 1.15 < r]^ < 2.44 with almost full 
azimuthal acceptance. The backward muon arm has a similar acceptance. We only 
consider leptons from the forward arm so as not to mix regions of 0:2. The ALICE ||23|| 
detector is designed to study nucleus- nucleus interactions at the LHC. They will mea- 
sure electrons at central rapidity, |?7e| < 0.9, and muons in a forward muon spectrometer 
with 2.5 <rj^<A. We impose the cuts described above on the lepton pseudorapidity 
in order to simulate detector acceptances. These are summarized in Table 0. However, 
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or o 






electron 0.0 < rje < 1.0 
muon 0.0 < 1]^^ < 1.0 


-0.35 < r/e < 0.35 
1.15 <ri^< 2.44 


-0.9 <r]e< 0.9 
2.5 <ri^< 4.0 



Table 2: The cuts on the pseudorapidity of single leptons for the three energy regimes 
studied are indicated both for electrons and muons. 

no cuts on minimum lepton energy or azimuthal angle were imposed. 

4 Results 

We begin by comparing the ratios of lepton pair cross sections with the input R^. 
Figure | shows the ratio of differential cross sections of e'^e'^ and /x^/i^ pairs from 
correlated DD and BB decays in pA and pp collisions at the SPS, RHIC and LHC 
energies as a function of the lepton pair invariant mass (solid curves). These results 

are compared to the ratios Rf{{x2), (Q)) (dashed curves) and R^{{x2), \J (Q^)) (dotted- 
dashed curves) at the average X2 and Q of each mass bin summed over y^j. All the 
results are integrated over the rapidity interval given in Table |. For brevity, we define 
the ratio of the cross sections as a function of lepton pair mass as 

_ daP^/dM 

Note that r„ follows i?^ to a very good approximation at all energies. The higher the 
energy, the better the agreement: at LHC r^r agrees with R^ very well. Note also that 

7Ai ' 



Rfi{x2), (Q)) ^ Rf{{x2), V(Q2)) in all cases. 

The results shown in Fig. ^ are obtained using a delta function for the fragmentation. 
If the Peterson function is used with = 0.15 and = {mc/rn}jYec [l^ , the outcome 



is very similar. The difference in due to the fragmentation function is smaller than 
or on the order of the statistical fluctuations in the Monte Carlo. 

The ratio always lies below R^ for two reasons. First, some fraction of QQ 
production is through the qq —>■ QQ channel. Quark shadowing effects are also included 
in the calculation oir„. Second, the phase space integration smears the shadowing effect 
relative to R^{{x2), {Q))- The qq contribution decreases with energy while Rf changes 
more slowly over the X2 region at the LHC, leading to better agreement at this energy. 

We have studied the influence of the qq channel by turning off the quark distribu- 
tions in our simulation, i.e. we flx = = for all quark flavours. The resulting 
ratio r„{f^ = 0) always lies above r^, indicating that quark effects reduce the cross 
section ratios. However, r„{f^ = 0) is still below R^, suggesting that other effects influ- 
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Figure 2: The ratio of e^e^ and js^^Ji^ pair cross sections from correlated DD and BB 
decays in and pp collisions, r^, defined in Eq. (pO|), at SPS, RHIC and LHC energies (solid 
curves). The nuclear gluon distribution, R^{x2,Q), at (^2) and (Q) (V (Q^)) of each M^j bin 
is indicated by the dashed (dotted-dashed) curves. 

ence r^-. For example, at the SPS with M = 1.5 GeV, ^ 1.12 and r^rif^ = 0) 1.16 
whereas i?^ ~ 1.19. 

In order to estimate the effect of the phase space integration, we reformulate the 
ratio of lepton pair cross sections as 



/ d{FS)QQRg{x2,Q)h{xi,X2,mT, . . .) 



(21) 

(i(PS)QQ/i(xi,X2,mr, • • •) 
where h{xi,X2,mT, ■ ■ ■) is a complicated function which contains the QQ production 
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cross section, fragmentation to mesons, and the subsequent meson decay to leptons. 
The integration over the QQ phase space is denoted by (i(PS)QQ. We define the average 
values 

(xa) = y c?(PS)QQX2/i(xi,X2,mT, ...), (Q) = / c?(PS)qq Q Xs, m^, • • •), (22) 

and expand as a Taylor series around (^2) and {Q) up to second order. The terms 
linear in {x2 — (3^2)) and {Q — (Q)) are zero after integration and the ratio becomes 

r„ = R {{X2},{Q)) + - "-^-^ {{X2-{X2))) 



2 dxi 

l d'Rf{{x2) ,{Q)) 
2 dQ^ 



dx2 dQ ^^^^ ~ (^v)iQ - (Qn) + ■■■ ■ (23) 

The sum of the correction terms is always negative, thus the ratio r^. always lies below 
Rf. A study of the partial derivatives (see Fig. 1) shows that d'^R^ /dx^ < in the 
vicinity of (^2) for each energy so that the first correction term is always negative. Near 
(g), we have d'^Rf/dQ^ > at SPS but < at RHIC and LHC. The first correction 
term dominates at SPS but the second one becomes dominant at the LHC energy. The 
net deviation of r„ from R^{{x2), (Q)) is then a sum of different terms refiecting both 
X- and Q-dependence of R^. As also seen in Fig. ||, the ratio deviates slightly more 
from R^{{x2), (Q)) for the electron pairs than for the muon pairs. The curvature of -R^ 
with X is stronger at larger values of x and due to the differences in rapidity coverage, 
the average values of X2 are larger for electron pairs. 

In order to estimate next to leading order (NLO) effects, we have also calculated 



the differential cross section of QQ production to NLO p4|. The results are shown in 
Fig. ^ where the ratio K = {da^^^ / dM) / {da^^ / dM) is given for pA and pp collisions. 
This K factor increases with M but decreases with energy. The figure also shows that 
the K factor for pA collisions, Ka, is approximately equal to that for pp collisions, K^^ 
since 

This result can be expected even though the K factor is large because the gg channel 
dominates both LO and NLO QQ production. The X2 regions probed at LO and NLO 
are also very similar. Equation (^) further implies that the ratio of the pA to pp cross 
sections at LO and NLO are approximately equal: 

(25) 
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Figure 3: The ratio of the NLO to LO differential cross sections for QQ production as a 
function of the invariant mass of the QQ pair. Left panel: cc pairs produced in p^ (solid) and 
pp (dashed) collisions at SPS (upper curves), RHIC (middle curves) and LHC (lower curves) 
energy. Right panel: The same for hb production at RHIC (upper curves) and LHC (lower 
curves). 



Thus, the NLO corrections to the ratios in Fig. || should be small. A similar result was 



found for Drell-Yan production at low masses ||25|, |26 



As explained earlier, the uncorrelated pairs might be used to obtain the ratio (-R^)^. 
If the nuclear effects turn out to be too small to be measured, that is, if ~ 1, the 
squared ratio might still be observable. We test whether the ratio of uncorrelated 
lepton pair cross sections, r^'^'^°" agrees with [R^Y the LHC where the probability 
of uncorrelated production is largest. Based on our previous discussion of the K factor, 
the NLO effects in the ratio r^'^™'''' are expected to remain small. 

The results for DD decays at the LHC are shown in Fig. It appears that the 
ratio of cross sections does follow the slope of (Q))]^. However, the difference 

between the ratios is slightly larger than in the correlated case. The origin of this 
difference can be understood by expanding two factors in the ratio of uncorrelated 
cross sections, r^^^^^ ^ as in Eq. ([23|) , and noting that the deviations (Q^) — (fff ^ind 
{xl^ — {x'lj'^ are larger for uncorrelated pairs. Integration over the rapidity of one of 
the leptons in a correlated pair increases the allowed phase space and thus the rapidity 
range of the remaining lepton. This greater available phase space directly increases the 
variation of X2 and Q for uncorrelated pairs, leading to the larger deviations between 
r™=- and [R^f. 

The slopes of the ratios r„ and R^ (see Fig. H) are determined by the momentum 
fractions probed by the interaction, as seen in a comparison with Fig. |l|. The average 
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Figure 4: Comparison of [R^{{x), (Q))]^ (dashed) and uncorrelated lepton pair cross sec- 
tions, r™'^™'^, for DD pair decays (solid) at the LHC. 



momentum fraction of the parton from nucleus, {X2), with muon pairs is 0.14 < {X2) < 
0.32 at the SPS, indicated in Fig. |]. In this x region, the ratio Rf{x, Q^) is decreasing. 
At RHIC, the range of momentum fractions probed is 0.003 < (^2) < 0.012, where 
the ratio is increasing quite rapidly. Finally, at the LHC, 3 x 10^^ < (^2) < 2 x 10""^ 
where Rf{x, Q^) is almost independent of x. The values of (^2) are typically larger for 
electron pairs at collider energies because the electron coverage is more central than 
the muon coverage, as seen in Fig. | and Table ||. 

We note here that the behaviour of Rf{x, Q^) depends not only on X2 but also on 
Q^. Therefore the slope of Rf{{x), (Q^)), seen in Fig. ||, does not directly follow any of 
the curves shown in Fig. [1| but evolves with from one curve to another. The average 
values (Q2) of cc decays are 7.58 < (Q^) < 48.5 GeV^ at the SPS, 9.46 < (Q^) < 141 
GeV^ at RHIC, and 11.4 < {Q^) < 577 GeV^ at the LHC. The corresponding values 
of {Q) are shown in Fig. |5|. The minimum average are larger for bb production, 
32.0 < (g2) < 54.3 GeV2 at RHIC and 37.9 < {Q^) < 156 GeV^ at LHC. 

The lepton pair cross sections from BB decays are also shown in Fig. ||. The results 
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Figure 5: We give (x2) and (Q) as a function of pair mass M for DD decays at the three 
energies. The soHd Hues are for fi^ fi^ pairs while the dashed is for e~e^ pairs. 

are quite similar to those from DD decays: at RHIC, the difference between r^j and 
is ~ 5% but decreases to 1 — 2% at the LHC. The slopes of these curves are different 
than for DD, they are almost constant with M. This constancy is due to the much 
narrower x and higher regions probed by hb production. In a narrow x region, the 
change in slope with x is small while at high Q^, evolution is slower, as seen in Fig. 
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Figure 6: Comparison of lepton pairs from correlated and uncorrelated DD and BB pairs 
at LHC (left) and RHIC (right). The differential number densities of e^e^ , e^/i^ and fJ-^j-i^ 
pairs are given as a function of M. The lepton pairs originate from correlated DD (solid), 
correlated BB (dashed), uncorrelated DD (dotted) and uncorrelated BB (dot-dashed) decays. 
The results are shown without K factors. 



Figure | shows the mass dependence of the differential number density distribution 



of lepton pairs in central collisions, obtained from Eqs. ([161) and ([TtI) through differ- 
entiation with respect to b and setting b = 0. These results are shown without any 
K factors, so that the final result will be shifted when a K factor is applied. In par- 
ticular, the relative normalization of the curves will change since the correlated pairs 
will be multiplied by K while the uncorrelated pairs are increased by K"^. Taking this 
correction into account, one sees that at the LHC the number of uncorrelated pairs at 
M > 2 — 3 GeV is of the same order as the correlated ones. At higher M, the uncor- 
related pairs are reduced and can be safely neglected, except in the e/i configuration 
where the uncorrelated pairs are larger. On the other hand, lepton pairs originating 
from BB decays are negligible at low M but become important at M > 5 — 6 GeV. At 
RHIC, leptons from correlated DD decays clearly dominate those from uncorrelated 
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pairs over all masses so that the uncorrelated pairs can be neglected. Note that at this 
energy also, leptons from BB decays are of the same magnitude as the correlated DD 
decays for M>4-5 GeV. 

In Ref. fl^ it was shown that, in AA collisions, lepton pairs from uncorrelated 
heavy quark decays tend to have higher masses than those from correlated decays. 
This is because the rapidity gap, Ay, between the quarks is larger when they are 
uncorrelated and oc 1 + cosh Ay. The uncorrelated rate is much smaller in pA 
interactions, however, which would reduce the probability of a large Ay relative to 
AA collisions. Including a finite rapidity cut severely reduces the number of leptons 
from uncorrelated decays falling within the acceptance. The effect of the rapidity 
cut is to reduce the average M for uncorrelated relative to correlated pairs. Thus in 
Fig. H lepton pairs from uncorrelated decays are most strongly reduced at large masses 
when the rapidity acceptance is small, as in the e"'"e~ coverage at RHIC. Since the 
e/i measurement involves mid-rapidity electrons and forward (or backward) muons, 
the rapidity coverage is largest and fewer uncorrelated decays are rejected. In fact, the 
larger e/i coverage at the LHC as well as the higher DD rates lead to more uncorrelated 
pairs at larger masses than correlated pairs. 

The results in Fig. || are for opposite-sign lepton pairs only. The uncorrelated yield 
can be determined using a like-sign subtraction since like-sign pairs are all trivially un- 
correlated. The number of uncorrelated opposite-sign pairs is equivalent to the number 
of like-sign pairs so it is relatively easy to estimate the real amount of uncorrelated 
pairs. After this estimate, the uncorrelated pairs can be used to study {Rf^, as shown 
in Fig. |. 

Finally we wish to mention that along with "normal" uncorrelated lepton pairs 
originating from DD and BB decays a number of "crossed events", i.e. uncorrelated 
leptons from DB and BD decays, are also produced. Their rate lies between those from 
uncorrelated DD and BB pairs. The total number of observed lepton pairs from heavy 
quark decays is the sum of all pairs from DD, BB, DB and DB decays. It is necessary 
to know the relative number of pairs produced in each case in order to decompose the 
results, as in Fig. ||. However, since uncorrelated BB pairs are generally negligible, 
we can also neglect these "crossed events" and assume that it is possible to deduce 
which lepton pairs came from correlated DD pairs. We do not consider any "crossed 
events" from D ov B decays combined with leptons from n or K decays. These should 
predominantly appear at lower masses and can be removed by like-sign subtraction. 



5 Conclusions 

We have calculated lepton pair production from heavy quark decays in pA and pp 
collisions using the EKS98 nuclear modifications to the parton distributions. We have 
shown that the resulting cross section ratios of pA to pp rates reflect the initial nu- 
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clear modifications of gluons rather well, as seen in Fig. ||. The deviations caused by 
phase space and the finite quark contributions are analyzed. Based on these results, 
we conclude that if the ratio of leptonic cross sections were measured within an ac- 
ceptable accuracy, the nuclear gluon distributions could be determined. Several such 
measurements have been proposed. 

The NA60 experiment ||2l| at the SPS plans to measure charm production through 
lepton decays. At this energy, a 20% antishadowing in manifests itself as a 10% 
enhancement in the ratio of cross sections due to the finite quark component. Assuming 
the measurements could be done to ~ 10% accuracy, NA60 results could be used to 
pinpoint the behaviour of nuclear gluon distribution in the range 0.17 ^ x ^OA. 

At higher energies, the gg channel becomes more important so that refiects the 
cross section ratio more closely. At RHIC, the PHENIX detector ||2^ will use lepton 
pairs to measure heavy quark production at x values as low as 3 x 10~^. In this x 
region, the dependence of gluon modifications should be larger than at the SPS. 
The nuclear effects vary from 20% shadowing for charm to slight antishadowing for 
bottom. The increased importance of the gluon channel reduces the difference between 
dCpA/df^pp and R^ to ~ 5%. 

The ALICE detector [2^ is designed to study nucleus-nucleus interactions at the 
LHC. The high LHC energy makes it ideal for probing nuclear effects of gluon distri- 
butions and studying saturation effects at low x. For x values as low as 3 x 10~^, the 
evolution in should be measurable, as shown in Fig. |l]. Here, the nuclear modifi- 
cations are an ~ 30% effect, and the difference between dapA/da^p and R^ is of the 
order of one percent. 

All three experiments, taken together, provide a large lever arm in {X2) and (Q^) 
over which to study the nuclear gluon distribution. Measurements at the same energy in 
p^ and pp interactions with high enough statistics should reveal much new information 
on R^. A better determination of R^ in turn provides valuable input to the study of 
AA interactions. 
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